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RESEARCH SUMMARY 


Potential gains from tree improvement of ponderosa pine in 
southern Idaho were assessed after 16 years of field testing. 
From a comparison of progenies of 268 half-sib families from 37 
populations growing at four test locations, genetic gains in 16- 
year height were related to (1) seed collections within seed zones, 
(2) seed collections from selected populations within seed zones, 
and (3) seed orchards developed after family trials. 


Analysis of population differentiation provided bases for re- 
evaluation of seed transfer guidelines. In general, seed for 
reforestation may be transferred from the collection area (1) 

+ 750 £t. (230 m)*elevation,«<(2) ;t.50 miles*(0.7 degrees) Latitude, 
and (3).+ 60 miles: (1.2 degrees) ‘longitude... /In addition; a igaim 
of about 4 percent in 16-year height can be expected from seed 
collections in populations of highest mean performance. 


Quantitative genetic analyses provided estimates of the 
additive genetic variance and phenotypic variance required for 
predicting gains in 16-year height after seed orchards are 
developed from family trials. Depending on the intensity of 
family selection and individual selection within families, genetic 
gains of 8 to 14 percent can be expected from seed orchards. 
Realization of these gains requires seed orchards to be composed 
of less than 150 genotypes. 
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INTRODUCTION 


The basis for a practical tree improvement program of ponderosa pine in southern Idaho is 
contained in 16-year old family evaluation tests. As detailed in preceding publications (Wang 
1967; Wang and Patee 1974, 1976; Rehfeldt 1980), primary objectives of the tests are: to 
develop seed zones and seed transfer guidelines that account for genetic adaptations to 
natural environments, to identify outstanding individual trees and populations from which 
wind-pollinated seeds will provide genetic gains to artificial reforestation, and to provide 
materials and information prerequisite to future breeding programs. 


The present analyses of 16-year height utilize experimental materials more fully described 
in a series of papers that have assessed performance at ages 3, 8, and 1l. Test plantations 
currently require thinning; therefore, it is appropriate to consider alternative plans for the 
future research program and for development of seed orchards. 


MATERIALS AND METHODS 


Wind-pollinated cones had been collected from 1 to 10 trees in each of 37 populations. 
Populations were separated by about 1 000 m in elevation and 2 degrees in latitude and longi- 
tude; they also represented habitat types that ranged from the relatively dry Pinus 
ponderosa/Purshia tridentata to the relatively moist Abies grandis/Vaccinium globulare. 
Thereby, much of the geographic distribution and ecological amplitude of the species in south- 
west Idaho was represented. Seedlings from 268 families were planted in four-tree row plots 
in each of 10 blocks at four planting sites. About 1.5 m separated trees within rows; 3 m 
separated rows. Test sites included Idaho City (1 200 m elevation), Holcomb (1 100 m), in the 
southern portion of the region; and Boulder Creek (1 450 m) and Jack's Creek (1 650 m) in the 
north. Since the last report of this series (Rehfeldt 1980), the plantation at Boulder Creek 
was thinned: the two shortest trees in each family plot were removed, and, at Jack's Creek, 
three replicates of exceptionally low survival were destroyed. 


Statistical analyses for assessing population differentiation and quantitative genetics 
were made on the 16-year height of individual trees. An appropriate analysis would involve 
calculation of main effects for sites and populations; nested effects of blocks in sites and 
families in populations; four interactive effects; and variance within plots. As explained 
previously, however, imbalances in the original design were accentuated by survival percentages 
that ranged from 22 to 80 percent at the four test sites. Because missing data were excessive, 
least squares analyses were required for accurate estimation of variance components. Yet, the 
large size of the matrix needed for a single analysis precluded inversion. Consequently, 
statistical analyses were completed in separate stages. 


Population differentiation was assessed from a least square analysis of variance (table 
1). One population, represented by a single family, was excluded from analyses. In addition 
to assessing differentiation, this analysis allowed calculation of least square means for 
populations within and among planting sites. These means were used in all subsequent compar- 
isons of population performance; they reflect the height expected if an equal number of 
observations had been available for each population in each block at all planting sites. 
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Table 1.--Form and results of analysis of variance of population performance for height at age 16 


Source of variance Degrees of freedom Components expected in each mean square Mean square Fi 
Sites 3 oF + 203.4305 1 + 1503.950% 489,559 SS Sm 
Blocks in sites 30 oF + 203.4308 1 144,829 189216 
Populations 35 op + 43.7405 6 + 174.2105 27,335 3.44** 
Populations x sites 105 of + 43.7405. 7,952 TS Sc 
Residual 6125 of 5,124 


**Statistical significance at the 1 percent level of probability. 


Patterns of population differentiation were related to physiognomic criteria of the seed 
source by multiple regression: 


Y.= bo +: byXy + boXo + b3X3 


Br 
where 
a least squares mean for population i on all sites as calculated in preceding analysis, 
bo = intercept, 
b, to b3 = regression coefficients, 
X, to X3 = independent variables of elevation (meters), latitude (degrees -40) and 


longitude (degrees -100), respectively. 


Effects of habitat types (recurring plant communities of potential climax) on population 
differentiation were also tested in multiple regression analyses. For these tests, the previous 
model was adjusted by including constant terms for each habitat type represented by the 36 popula- 
tions. But, as shown previously (Rehfeldt 1980), inclusion of the effects of habitat types did 
not improve the fit of the data to the model. Habitat types apparently have little effect on 
differentiation and are subsequently ignored. 


The regression model presented above was also used on least squares mean height of popula- 
tions from each of the four sites. As detailed later, however, the mean height of populations 
between ages 11 and 16 was so strongly correlated among sites that the results of the separate 
regression analyses essentially duplicated those presented previously (Rehfeldt 1980). Conse- 
quently, results of separate analyses for each site are not presented. 


Guidelines for seed transfer were developed from patterns of population differentiation. 
According to techniques used previously (Rehfeldt 1979), the analysis of variance allows 
calculation of the least significant difference (lsd) among population means that provide for 
statistical detection of mean differences at given levels of probability (Steele and Torrie 
1960). A value of lsd was calculated at the relatively low level of probability of 0.2 in 
order to guard against making the error of accepting no differences among populations when 
differences actually exist. Then, the expression Isd(9.5)/b, where ‘b is the regression 
coefficient, provides the distance in elevation, latitude, and longitude associated with mean 
differences at the 80 percent level of probability. These distances are proposed as guidelines 
for limiting seed transfer. 


Possible genetic gains within seed zones were estimated from quantitative genetic analyses 
made for each site according to a model of random effects (table 2). These analyses allowed 
estimates of family heritabilities, individual heritabilities, and genetic gains expected 
after alternative plans of seed orchard development. Estimates of genetic gains from family 
selection and individual selection within families followed Namkoong (1979). 
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Table 2.--Form of analyses of variance for quantitative genetics 


Source of variance Expected mean squares? 
Blocks OW ko 15 + k_ 02 
1) 5.9uB 
o2 2 > 2 
Populations Wat ky Geert ke CF /p + ky op 
eat A ° o2 vy) 2 
Families in populations W + ky OF ky OF /p 
os y 
Experimental error! W+ ky oF 
Within ow 
‘Contains all sources of variance involving interaction of blocks 
2where: 
Holcomb Idaho City Boulder Creek Jack's Creek 
k, = 1.88 1.90 1.89 1556 
k= 11.36 16.04 16.29 4.70 
kes 1201 16.83 16.75 Dai 
Kis 74.50 109.04 110.47 30.45 
n= 334.56 377.89 397.03 177.01 


RESULTS AND DISCUSSION 
Population Differentiation 


Analyses of population differentiation readily detected mean differences among populations 
across,all planting sites (table 1). But, main,effects.of.sites.were. also significant, and 
populations tended to perform differently at the various planting sites. The interaction 
apparently arises from expression of differential adaptations of populations at planting sites 
that encompass a large elevational (550 m) and latitudinal (> 1 degree) range. 


Least square mean heights and physiognomic criteria of the seed sources are presented in 
table 3. These means are representative of mean performance across all four test sites, and, 
therefore, they are confounded by interactions between genotype and environment. But, simple 
correlation coefficients (r) of mean height for populations at each site and least square 
means ranged from 0.64 to 0.85. The high values of these correlation coefficients suggest 
that confounding of means by genotype-environment interactions is small. Least square means 
apparently reflect an average performance suitable for assessing differentiation in relation 
to physiognomic criteria of the seed source. But, because the correlations were not perfect, 
differentiation also was assessed from performance at each site. 


Table 3.--Location and least square mean height of populations according to origin 


Popula- 
tion Location Latitude Longitude Elevation 16-year height 
Number Degrees Degrees Meters cm 
1 Anderson 44.55 116.20 1520 292 
2 Mill Creek 44.45 116.30 1370 294 
5 Camp Creek 44.62 115.65 1830 286 
6 Crawford 44.53 1155.97 1460 296 
9 Warm Springs 
Ridge 43.85 115.88 1520 296 
iat Lester Creek 43.62 115.38 1830 , 281 
14 Sagehen 44.32 11G;:2'3 1580 308 
15 Third Fork 44.37 116.33 1250 309 
16 Powellson Creek 44.78 115.87 1770 293 
17 Circle Creek 45.05 116.527 1280 311 
18 Bear Creek 43.95 115.47 1830 » 275 
19 Idaho City 43.672 UTS. 72 1460 301 
22 South Fork 
Payette River 44.07 115.50 1220 299 
23 Mill Creek 44.67 116.87 1250 309 
24 Pine Creek 44.70 116.83 1160 309 
25 Babbit Creek 43.80 LS 272 1950 283 
26 Mica Creek 44.60 116627 1520 285 
27 Cabin Creek 44.63 116.27 1370 331 
28 Mud Creek 45.05 116.38 1520 308 
29 Rough Creek 44.90 116.48 975 314 
35 Dutch Creek A 43.80 115.37 1520 289 
38 Gircle- "Cc" 45.05 116.27 1220 318 
39 Meadow Valley 44.98 116.22 1220 318 
41 Dutch Creek B 43.73 315.45 1980 288 
42 Trail Creek 43.60 15.275 1190 315 
46 West Fork Creek 44.30 11587 1370 286 
47 Garden Valley 44.07 115.93 1040 315 
50 Bumgartner 43.62 115.08 1310 271 
5S Camp Creek 44.88 115.70 1220 320 
54 Zena Creek 45.07 115.75 1520 271 
20A Idaho City 43.77 11S..72 1280 296 
20B Idaho City 43.87 115.583 1280 282 
21A Hazard Creek 45.18 116.30 1070 321 
21B Scriver Creek 44.17 116.00 1310 318 
43A Rice Creek 45.75. 1163.33 1430 301 
43B Fleming Creek 44.03 116.10 1310 296 


In the multiple regression analysis, elevation, latitude and longitude of the seed source 
accounted for 51 percent of the variation in mean performance of populations (table 4). 
Visual examination of the residuals from regression suggested that transformations of the 
independent variables would not necessarily improve the fit. The regression statistics show 
that for an elevational gradient of 1 000 m, mean height of populations decreases by about 10 
percent; and across a geographic interval of one degree of latitude or longitude, mean height 
changes by 3 percent and 2 percent, respectively. 


Since the simple correlations (r) of population means at age 16 with those at age 11 were 
0.78, 0.89, 0.93, and 0.91 at Jack's Creek, Holcomb, Boulder Creek, and Idaho City, respectively, 
it is not surprising that multiple regression models that fit mean performance at each site to 
physiognomic variables essentially duplicated results of similar models for height at age 11 
(Rehfeldt 1980). Moreover, models that involved mean values of 16-year height across all 
sites (table 4) accounted for more variance than those models involving 1l- or 16-year height 
at a single site. Undoubtedly, analyses that combined observations from all sites provided 


better estimates of population means and greater statistical power than previous analyses that 
were made for each site separately. 


Table 4.--Results of multiple regression analyses for relating population performance 
to physiognomic criteria of the seed source 


Variable Regression coefficient Isd(9 2) 
(b) b 

Elevation (X)) -0.026 466 m (1500 ft) 
Latitude (X>) 8.791 1.4° (100 miles) 
Longitude (x3) 5.134 2.3° (120 miles) 

by = 268.84 

R2 = 0.51 

eat & =) 10899 

yx 


Statistical analyses have detected differences among populations in 16-year height and 
have related differentiation to geographic and ecologic conditions of the seed source. The 
applicability of these results to practices of seed zoning depends on the degree by which tree’ 
height reflects adaptations to the natural environment. While there is little doubt that seed 
zones should be based on numerous adaptational characteristics, no data other than height are 
available for ponderosa pine in southern Idaho. But, adaptations in ponderosa pine may be 
manifest much like adaptations in Douglas-fir of north Idaho: the growth, phenology, and cold 
hardiness of populations are so highly intercorrelated that variation among populations in a 
single character reflects numerous characters (Rehfeldt 1979). Consequently, it is assumed 
that observed patterns of variation in 16-year height of ponderosa pine are adaptive. 


Table 4 reveals the distance in elevation, latitude, and longitude that is associated 
with mean differences among populations equal to a value of lsd(g.2). At minimum distances of 
1,533 ft (466 m) elevation, 100 miles (1.4 degrees) latitude, or 120 miles (1.9 degrees) 
longitude, differences among populations are detectable with a probability of 80 percent. Or 
as a corollary, seed for afforestation may be transferred from the collection area: + 750 ft 
(230 m) elevation, + 50 miles (0.7 degrees) latitude and + 60 miles (1.2 degrees) longitude. 


Yet, it should be recognized that the recommended seed zones are based solely on the 
probability of detecting differences in height among populations. Not only do these zones 
fail to account for traits other than height, they also fail to account for administrative 
factors such as patterns of land ownership or location of future planting programs. However, 
these guidelines can be altered readily provided that an increased risk of maladaptation is 
recognized. A difference in height of about 10 percent can be expected from populations 
separated by 1 000 m (table 4); populations separated by 1 degree of latitude or longitude can 
be expected to,difter by..35. percent and 2 percent, respectively. 


Genetic gains from tree improvement are readily attained by collecting seeds for refores- 
tation from natural stands of proven genetic value. For instance, if the present region is 
divided into four zones separated at 1 400 m elevation and 44.3 degrees latitude, a gain in 
height of about 4 percent could be expected within zones when seed collections are made from 
the following populations: 


Zone Superior populations* 
Low north 2739.5 55)20A 
High north LAK 28 
Low south 42,47,21B 
High south 14** 19 


*Coded in table 3. 
**Suitable for both northern and southern zones. 
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Quantitative Genetics 


All mean squares associated with the main effects of families at each planting site 
(table 5) were statistically significant at the 1 percent level of probability. As a percent 
of the sum of variance components, the family variance accounted for about 4 percent of the 
total variance and slightly more variance than effects of populations. Still, the use of 
these results to estimate genetic components of variance (table 6) needs explanation. 


Because of the methods of calculation, bias is introduced at four levels into estimates 
of the additive genetic variance, phenotypic variance, and, consequently, heritabilities. 
First, the assumption is implied that families were composed of half-sibs. But, progenies of 
wind-pollinated trees are slightly more closely related than actual half-sibs (Namkoong 1966); 
thus, additive genetic variances (table 6) are slightly overestimated. 


Table 5.--Results of analyses of variance of 16-year height for each location 


Holcomb Idaho City Boulder Creek Jack's Creek 

Source of Mean Mean Mean Mean 
variance d.f. square o2 d.f. square Ge d.f. square o*  d.f. square o2 
Blocks 7 544,147 1604 9 213,488 1037 9 154,743 363 5 178,051 945 
Populations 36 19,786 127 36 21,5302 OF, 36 30,265 128 36 22,499 315 
Families/ d 

Populations 198 10,350 249 207 10,689 218 207 16: 171 340 193 12,907 464 
Experimental 

error! 1186 7,927 2002 1826 7,198 1924 1852 10,632 3804 470 10,726 2338 
Within 2041 3.5762 3762 2948 3,542 3542 2095 3,442 3442 710 7,078 7078 


IContains sources of variance involving interactions of blocks. 


Table 6.--Genetic components of variance for each location. 
Calculations follow Namkoong (1979) 


Holcomb Idaho City Boulder Creek Jack's Creek 


Additive genetic variance: 


oe 2 
on > 405 7p 996 872 1360 1856 


Total phenotypic variance: 


+ 02 6013 5684 7586 9880 


aes 2 
=o +o E/P 


Family phenotypic variance: 


992 661 1000 2746 


F/P 


Bese se AG 0.27 0.33 0.34 0.17 


Heritability of individuals: 
2 

en 
of 
eo EEE Ee 
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Secondly, since estimates of additive genetic variance from single planting sites are 
confounded by interactions of genotype and environment (Namkoong and others 1966), present 
estimates (table 6) may be subject to bias. But, genetic gains for ponderosa pine in southern 
Idaho will accrue within seed zones, and, therefore, appropriate estimates of heritabilities 
should be inflated by the effects of the genotype-environment interaction (Namkoong 1979). 


Thirdly, estimates of the phenotypic variance apparently do not contain components of 
variance due to genotype-environment interaction. It is likely, however, that some of the 
effects of a genotype-environment interaction expected within seed zones are contained within 
the interaction of blocks and families. Whereas the mean height of trees varied 63 cm among 
planting sites, mean differences associated with blocks ranged from 67 cm to 105 cm at the 
various planting sites. Significant mean differences among blocks (table 5) arise partially 
because blocks were arranged according to aspect at each site. Thereby, blocks reflect a 
variety of sites typical of the seed zone represented by each planting site. Therefore, the 
interaction among blocks and families (table 5) likely reflects genotype - environment inter- 
actions within seed zones as well as sampling errors. 


Finally, sampling errors are magnified by the experimental design. Only four trees were 
planted in each block at each test site. As reflected by harmonic means (k, in table 1), an 
average of two or fewer trees remain. Thus, family plots poorly reflect family means in each 
block; experimental errors (block interactions) are magnified; and phenotypic variances are 
greatly exaggerated in relation to the additive genetic variance. 


Thus, there is little doubt that the genetic components of variance (table 6) are only 
approximate. But, particularly if each planting site is considered as a representative of a 
single seed zone, errors of estimate can be absorbed without invalidation of statistical 
analyses. 


Results of analyses of variance (table 5) and calculations of genetic components of 
variance (table 6) are similar to those from age 11 (Rehfeldt 1980). Most of this corre- 
spondence can be traced. to the high correlation of individual tree heights (r =.0.86, 0.88, 
0.91, and 0.83 at Jack's Creek, Boulder Creek, Idaho City, and Holcomb, respectively) and 
family means (r = 0.81, 0.88, 0.89 and 0.81, respectively) between data at age 11 and age 16. 
In fact, analyses of variance of the deviation from regression of 16-year on ll-year height 
detected no differences among blocks, families, or their interaction; all of the variance was 
attributable to experimental errors or to variation within family plots. Nevertheless, as 
observed for ponderosa pine in California (Namkoong and Conkle 1976), the proportion of the 
total variance attributable to families within populations has decreased between ages 8, 11, 
and 16. This reduction occurred because error variances (within plots and experimental 
error) increased. In fact, while error components of variance doubled between ages 11 and 16, 
family variance components increased in absolute size by only 50 percent. Consequently, 
estimates of additive genetic variance are less; phenotypic variances are greater; herita- 
bilities are less; and expected gains are smaller at age 16 than at age 11. 


Genetic gains in 16-year height that are expected after interpollination of trees chosen 
under various intensities of family and individual tree selection are detailed in figure 1. 
These gains were calculated for family selection by: 


Pets i) 
AG ihe Onp 


and for selection of individuals within selected families by: 


E wee 
AG 0.75ih¢ On 


where 


I 


AG = genetic gain in units of measure 


i = selection intensity; a constant relation to the proportion of the population 
selected (Namkoong 1979). 


Zia . ‘ 
hes hy, op and Opp are defined in table 6. 
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Figure 1.--Expected gains (ordinate) 


0 0 in 16-year height from alternative 
337-20. 10 33 20 10 levels of family selection (bars) and 
individual selection (shading) with- 
PERCENT FAMILY SELECTION in families for*each’ planting srte: 


Gains presented in figure 1 are expressed as a percentage of the mean height of only 
those trees at each planting site that represent the seed zone within which that site lies. A 
3 to 5 percent gain over local seed sources can be expected after intermatings of the tallest 
families. And, gains of 10 to 15 percent seem likely in progenies of the best performing 
individual trees within selected families. Relatively high gains will be realized, however, 
only if seed orchards are composed of relatively few genotypes (fig. 2). And therefore, an 
additional investment will be required to develop clonal seed orchards. 


These results attest to potential genetic gains at several levels of tree improvement for 
ponderosa pine. First, seed transfer guidelines account for natural adaptation of populations 
and reduce losses in productivity from maladaptation. These guidelines limit seed transfer to 
+ 750 £t (230 m) elevation, + 50 miles (0.7 degrees) latitude, and + 60 miles (1.2 degrees) 
longitude. Secondly, a 4 percent gain in 16-year height can be attained by collecting seeds 
for reforestation from populations of proven genetic potential. And, thirdly, depending on 
the intensity of selection, gains from 8 to 14 percent can be expected in the next generation 
from seed orchards developed from select individuals within selected families. 
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Figure 2.--Approximate relationship between per- 
cent gain in 16-year height and number of geno- 
types selected for seed orchards at each plant- 
ing site. The number of genotypes was derived 
from the various combinations of family and 
individual tree selection detailed in figure 1. 
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The 16-year heights of ponderosa pine progenies from 268 
half-sib families representing 37 southern Idaho populations 
were compared on four test sites. Analyses of population 
differentiation showed that a genetic gain in 16-year height 
of about 4 percent will result from collections of seed in 
selected populations. But seed for reforestation should not 
be transferred more than 460 m elevation, 1.4° latitude or 
2.4° longitude. Quantitative genetic analyses implied that 
genetic gains of 8 to 14 percent can be expected in the height 
of 16-year-old progeny of selected families and individual 
trees. 
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The Intermountain Station, headquartered in Ogden, 
Utah, is one of eight regional experiment stations charged 
with providing scientific knowledge to help resource 
managers meet human needs and protect forest and range 
ecosystems. 

The Intermountain Station includes the States of 
Montana, Idaho, Utah, Nevada, and western Wyoming. 
About 231 million acres, or 85 percent, of the land area in the 
Station territory are classified as forest and rangeland. These 
lands include grasslands, deserts, shrublands, alpine areas, 
and well-stocked forests. They supply fiber for forest in- 
dustries; minerals for energy and industrial development; and 
water for domestic and industrial consumption. They also = 
provide recreation opportunities for millions of visitors each 
year. 


Field programs and research work units of the Station 
are maintained in: 
Boise, Idaho 


Bozeman, Montana (in cooperation with Montana 
State University) 


Logan, Utah (in cooperation with Utah State 
University) 


Missoula, Montana (in cooperation with the 
University of Montana) 


Moscow, Idaho (in cooperation with the Univer- 
sity of Idaho) 


Provo, Utah (in cooperation with Brigham Young 
University) 


Reno, Nevada (in cooperation with the University 
of Nevada) 


ni} . } 
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